What is claimed is: 

1. A method for detemiining candidates for gene deletions and additions using a 
model of a metabolic network associated with an organism, the model comprising a 
plurality of metabolic reactions defining metabolite relationships, the method comprising: 

5 selecting a bioengineering objective for the organism; 
selecting at least one cellular objective; 

forming an optimization problem that couples the at least one cellular objective with the 

bioengineering objective; and 
solving the optimization problem to yield at least one candidate. 
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2. The method of claim 1 further comprising modifying the organism with the 
candidate. 

3. The method of claim 1 wherein the bioengineering objective is overproduction of a 
15 chemical. 

4. The method of claim 1 wherein the bioengineering objective is underproduction of 
a chemical. 

20 5. The method of claim 1 wherein the cellular objective is growth. 

6. The method of claim 1 wherein the cellular objective is minimization of metabohc 
adjustment. 

25 7. The method of claim 1 wherein the candidate is a candidate for gene deletion, and 
the optimization problem includes a binary value for specifying if reaction flux is active or 
inactive. 

8. The method of claim 1 wherein the optimization problem is a bilevel optimization 
30 problem. 
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9. The method of claim 1 wherein the optimization problem is a mixed-integer 
optimization problem. 

10. The method of claim 1 wherein the optimization problem includes at least one 
stoichiometric constraint. 

11. The method of claim 1 wherein the optimization problem includes at least one 
chemical uptake constraint. 

12. The method of claim 1 wherein the step of forming an optimization problem 
includes quantifying the cellular objective as an aggregate reaction flux. 

13. The method of claim 1 further comprising evaluating performance limits of the 
metabolic network with the at least one candidate based on ability of the network to meet 
the at least one cellular objective. 

14. The method of claim 1 wherein the cellular objective is selected from the group 
consisting of: maximizing a growth rate, maximizing ATP production, minimizing 
metabolic adjustment, minimizing nutrient uptake, minimizing redox production, 
minimizing a Euclidean norm, and combinations thereof. 

15. The method of claim 1 wherein the bioengineering objective is overproduction of 
glycerol and at least one candidate is for gene deletion and comprising genes coding for the 
enzymes fructose- 1,6-bisphosphatase, fructose- 1,6-bisphosphatase aldolase, 
phosphoglycerate kinase, glyceraldehydes-3-phosphate dehydrogenase, 
phosphoenolpyruvate synthase, NADH dehydrogenase I, phosphogluconate dehydratase, 2- 
keto-3-deoxy-6-phosphofluconate aldolase, triosphosphate isomerase, glucose 6- 
phosphate-1 -dehydrogenase, 6-phosphogluconolactonase, deoxyribose-phosphate aldolase, 
aldehyde dehydrogenase, or combinations thereof. 
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16, The method of claim 1 wherein the bioengineering objective is overproduction of 
1,3-propanediol and at least one candidate is for gene deletion and comprising genes 
coding for the enzymes fructose- 1,6-bisphosphatase, fructose- 1,6-bisphosphatase aldolase, 
phosphoglycerate kinase, glyceraldehyde-3-phosphate dehyrogenase, triosphosphate 

5 isomerase, glucose 6-phosphate-l -dehydrogenase, 6-phosphogluconolactonase, 
deoxyribose-phosphate aldolase, aldehyde dehydrogenase, or combinations thereof. 

17. The method of claim 1 wherein the bioengineering objective is overproduction of 
succinate and at least one candidate is for gene deletion and comprising genes coding for 

10 the enzymes pyruvate formate lyase, acetaldehyde dehydrogenase, pyruvate kinase, FOFl- 
ATPase, NADH dehydrogenase I, fumarase, D-Lactate dehydrogenase, pyridine nucleotide 
transhydrogenase, phosphotransacetylase, acetate kinase, phosphotransferase, or 
combinations thereof. 

15 18. The method of claim 1 wherein the bioengineering objective is overproduction of 
lactate and at least one candidate is for gene deletion and comprising genes coding for the 
enzymes phosphotransacetylase, acetate kinase, phosphofructokinase, fructose- 1,6- 
bisphosphatase aldolase, triosphosphate isomerase, acetaldehyde dehyrogenase, 
glucokinase, or combinations thereof. 
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19. A computer-based method for determining candidates for gene deletions and 
additions using a model of a metabolic network associated with an organism, the model 
comprising a plurality of metabolic reactions defining metabolite relationships, the method 
comprising: 

25 inputting at least one bioengineering objective; 
receiving as input as least one cellular objective; 

forming an optimization that quantifies the cellular objective as an aggregate reaction flux 

and couples the at least one cellular objective with the bioengineering objective; 
solving the optimization problem to yield at least one candidate; and 
30 outputting the at least one candidate. 
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